The quantum walk temperature 
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A thermodynamic theory is developed to describe the behavior of the entanglement between the 
coin and position degrees of freedom of the quantum walk on the line. This theory shows that, 
in spite of the unitary evolution, a steady state is established after a Markovian transient stage. 
This study suggests that if a quantum dynamics is developed in a composite Hilbert space (i.e. the 
tensor product of several sub-spaces) then the behavior of an operator that only belongs to one of 
the sub-spaces may camouflage the unitary character of the global evolution. 
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I. INTRODUCTION 

The concept of isolated system plays a fundamental 
role in the formulation of the quantum mechanics. This 
concept is an idealization that was constructed as an aid 
to understand some phenomena displayed by real sys- 
tems which may be regarded as approximately isolated. 
However, since about 50 years ago, the study of quan- 
tum decoherence has acquired a central position in the 
formulation of the quantum mechanics. In fact, concepts 
such as thermodynamic equilibrium seem impossible to 
coordinate with the idea of isolated system because the 
quantum state for such a system follows a unitary evo- 
lution and it cannot reach a final equilibrium state at 
t — > oo. 

In this context we ask ourselves if it is possible to intro- 
duce the concept of temperature for an isolated quantum 
system which evolves in a composite Hilbert space. In 
this paper the system known as the quantum walk on 
the line [l[ has been chosen as a model to answer this 
question . The quantum walk (QW) is a natural gen- 
eralization of the classical random walk in the frame of 
quantum computation and quantum information process- 
ing and it is receiving permanent attention 043- It has 
the property to spread over the line linearly in time as 
characterized by the standard deviation a(t) ~ t, while 
its classical analog spreads out as a(t) ~ t 1 / 2 . This 
property, as well as quantum parallelism and quantum 
entanglement, could be used to increase the efficiency of 
quantum algorithms [H, 0] . Recently we have been inves- 
tigating 0-Q the asymptotic behavior of the QW on the 
line, focusing on the probability distribution of chirality 
independently of position. We showed that this distribu- 
tion has a stationary long-time limit that depends on the 
initial conditions. This result is unexpected in the con- 
text of the unitary evolution of the QW because such a 
behavior is usually associated to a Markovian process. In 
this paper we further explore the behavior of the chiral- 



ity distribution and define a thermodynamic equilibrium 
between the degrees of freedom of position and chirality. 
This equilibrium allows to introduce a temperature con- 
cept for this unitary closed system. We obtain a master 
equation with a time-dependent population rate, that de- 
scribes the transient behavior of the reduced density op- 
erator of the QW towards thermodynamic equilibrium. 
The QWs reduced density operator shows an surprising 
behavior: Its behavior looks diffusive but however the 
global evolution of the system is unitary. 

The paper is organized as follows. In the next section 
the standard QW model is developed, in the third sec- 
tion the entanglement temperature is defined. Then in 
the fourth and fifth sections the entanglement temper- 
ature is obtained for localized initial conditions and for 
distributed initial conditions respectively. In the six sec- 
tion the transient behavior towards thermal equilibrium 
is studied, and finally in the last section some conclusions 
are drafted. 



II. QW ON THE LINE 

The composite Hilbert space of the QW is the tensor 
product Hs^Hc where V. s is the Hilbert space associated 
to the motion on the line and T-L c is the chirality (or coin) 
Hilbert space. In this composite space the walker moves, 
at discrete time steps t € N, along a one-dimensional 
lattice of sites k e Z. The direction of motion depends 
on the state of the chirality, with the eigenstates R and 
L. The wave vector can be expressed as the spinor 
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where the upper (lower) component is associated to the 
left (right) chirality. 

Then P kL {t) = \a k {t)\ 2 and P kR (t) = \b k {t)\ 2 denote 
the probability of finding the walker at (fc, t) and the 
coin in state R and L, respectively. The probability of 
finding the walker at (k, t) is 
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and£ fe P(M) = l- 

The QW is ruled by a unitary map whose standard 
form is [10h13| 

a k (t +l) = a k+1 (t) cos 9 + b k+1 (t) sin 9, 

b k (t + 1) = a k -i(t) sine - b k -x(t) cosO. (3) 

where 9 € [0, tt/2] is a parameter defining the bias of the 
coin toss (9 = j for an unbiased or Hadamard coin). The 
global left and right chirality probabilities is defined as 



This interesting result for the QW shows that the long- 
time probability to find the system with left or right chi- 
rality has a limit. Therefore, although the dynamical 
evolution of the QW is unitary, the evolution of its GCD 
has an asymptotic limit characteristic of a diffusive be- 
havior. This situation is further surprising if we compare 
our case with the case of the QW on finite graphs [l4j 
where it is shown that there is no converge to any sta- 
tionary distribution. 
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with Pji(t) +Pt(t) = 1. The global chirality distribution 
(GCD) is defined as the distribution formed by the couple 

pkj^x . It is shown in Ref. [7] that the GCD satisfies 

the following map 
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where 



Q(t)= J2 a k{t)b* k {t). 



(5) 



(6) 



The two dimensional matrix in Eq. ([5} can be interpreted 
as a transition probability matrix for a classical two di- 
mensional random walk as it satisfies the necessary re- 
quirements, namely, all its elements are positive and the 
sum over the elements of any column or row is equal to 
one. On the other hand, it is clear that Q(t) accounts for 
the interferences. When Q(i) vanishes the behavior of the 
GCD can be described as a classical Markovian process. 
However in the generic case Q(t) together with Pt(i) 
and P R (t) are time depend functions that have long-time 
limiting values @ which are determined by the initial 
conditions of Eq.©. Eq.© can be solved in this limit. 
We define 
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lim P L (t) 
t — > oo 

lim P R {t) 
t — > oo 

lim Q(t) 
t -> oo ' 



(7) 



and then we obtain the asymptotic stationary solution 
for the GCD as 
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1 + 2Re(Q )/tan( 
1 - 2Re(Q )/tan( 



(8) 



III. ENTANGLEMENT AND TEMPERATURE 

The concept of entanglement is an important element 
in the development of quantum communication, quantum 
cryptography and quantum computation. In this context 
several authors have investigated the relation between 
asymptotic entanglement and the initial conditions of the 
QW [151 - 4251 ] . Other authors [26|, |27J have proposed to use 
the QW as a tool for quantum algorithm development 
and as an entanglement generator potentially useful to 
test quantum hardware. 

The unitary evolution of the QW generates entan- 
glement between the coin and position degrees of free- 
dom. To characterize this entanglement we start with 
the von Neumann entropy which is quantum analogue of 
the Gibbs entropy 



Sn{p) = -tr(plogp). 



(9) 



where p — \^ (t)) (t)\ is the density matrix of the quan- 
tum system. Due to the unitary dynamics of the QW the 
system remains in a pure state and this entropy vanishes. 
However for these pure states the entanglement between 
the chirality and the position can be quantified by the 
associated von Neumann entropy for the reduced density 
operator [UlUll that defines the entropy of entanglement 



where 



S(p) = -trf> c logp c ), 



Pc = tr(p), 



(10) 



(11) 



and the partial trace is taken over the positions. Using 
the wave function Eq.© and its normalization properties 
we obtain the reduced density operator 15] 
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The density operator p c has the eigenvalues 



A+ = - 



1 ± 



^l-4(p L (t)P R (t)-\Q(t)f 



(12) 
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Then the entanglement entropy Eg. (1111) is expressed 
through these two eigenvalues as 



S(P) 



-A+ log A + — A_ log A_ 



(14) 
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In the asymptotic regime t — > oo the eigenvalues go to 
stationary limit, X± — > A± and 



A+ = 



l±i/l- 



4(n L Ur-\q \ 2 ) 



Further from Eq.© follows the relation 



(15) 



(16) 



which is substituted in Eq. (115[) . and then the asymptotic 
eigenvalues are expressed as 



with 



A± = - ± VX. 



x = |Qo| 2 + (Re(Q o )/tan0) 2 . 



(17) 



(18) 



Note that the values of the interference term Qo are con- 
strained to satisfy the condition 



< A+A_ < 1, 



and then 



< X < 
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FIG. 1: Dimensionless thermodynamic function normalized 
by log(2) as a function of the dimensionless parameter \- 
From the top to the bottom they are: in thick line f3e, in 
dashed line the entropy So, in thin line the energy /3U and in 
dashed-dot line the Helmholtz free energy ft A. 



The entanglement entropy has an asymptotic limit too 

S a = -A+logA + - A_logA_, (21) 

that only depends on the initial conditions through the 
interference term Qq. Therefore we are led to consider 
that after some transient time the QW achieves a thermo- 
dynamic equilibrium between the position and chirality 
degrees of freedom. In order to make a fuller description 
of this equilibrium it is necessary to connect the eigenval- 
ues of p c with its associated Hamiltonian operator H c . To 
obtain this connection we shall use the quantum Brow- 
nian motion model of Ref. 28[. We considered the sys- 
tem associated with the chirality degrees of freedom and 
characterized by the density matrix p c in thermal contact 
(entanglement) with the bath system associated with the 
position degrees of freedom, the lattice. In this context 
p c satisfies the equation 



^£ - —m 
at ~ m l c 



(22) 



where [H c ,p c ] is the commutator and Tp c represents the 
Brownian motion of p c induced by the noise (fluctuating 
forces) exerted on p c by the lattice (position degrees of 
freedom). In the equilibrium (stationary) situation we 
must have dp c /dt = and Tp c = [Hj], that is 



\H r . 



= 0. 



(23) 



Therefore, in the asymptotic regime, the density opera- 
tor p c must be an explicit function of the Hamiltonian 



operator which must be time independent. Now we call 
{$+,$_} the eigenfunctions of the density matrix, and 
then the operators H c and p c are both diagonal in this ba- 
sis. Therefore, the eigenvalues A + and A^ depend of the 
corresponding eigenvalues of H c - We take these eigen- 
values to be {— e, e} without any loss of generality, they 
represent the possible values of the entanglement energy. 
This interpretation agrees with the fact that A + and A_ 
are the probabilities that the system is in the eigenstate 
<I> + or The precise dependence between A± and ±e 
is determined by the type of ensemble we construct. The 
main proposal of this paper is that this equilibrium cor- 
responds to a quantum canonical ensemble. Therefore 
we propose that 



A+ = 



,±/9e 



-/3e> 



(24) 



which defines the entanglement temperature T = 1//3. 
Then the probability that the state chosen at random 
from the ensemble {$+, possesses an energy e is 

determined by the Boltzmann factor e~ l3e . Let us call p c 
the diagonal expression of the density operator p c , then 



A+ 
A_ 



1 



;/3£ _|_ e -f3e 



e /3e 






-Pe 



(25) 



This operator is formally the same density operator that 
corresponds to an electron with possesses an intrinsic spin 
and a magnetic moment in a external magnetic field [31] . 
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Starting from Eq. ([25|) it is possible to build the ther- 
modynamics for the QW entanglement. The partition 
function of the system is then given by 



2cosh(/3e). 



(26) 



Accordingly, and also using Eqs. (|151l2"4l . the temperature 
is given by 



T = 2e/ln 



1 



i- VxV 



the Helmholtz free energy by 
.4 



iln[2cosh(/3e)] = |lnQ 



the internal energy by 

U = -etanh(/3e) = -2e^/x, 
and finally the entropy by 

S = PU- BA, 



(27) 



(28) 



(29) 



(30) 



where this last thermodynamic definition for the entropy 
of course agrees with the previous Shannon expression in 
Eq.(f2~Tj). To finished this section in Fig. [T]we present the 
dependence of these thermodynamic magnitudes with the 
interference parameter \- 



IV. LOCALIZED INITIAL CONDITIONS 

As seen in the previous section the thermodynamics of 
entanglement only depends on the interference term Qo 
which in turn only depends on the initial conditions, as 
shown in Q- 

In order to investigate this dependence on the initial 
conditions of the system we consider first the localized 
case. The initial state of the walker is assumed to be 
sharply localized at the origin with arbitrary chirality, 
thus 



l*(0)> = 



cos(7/2) 
expitp sin(7/2) 



|0), 



(31) 



where 7 6 [0, tt] and <p € [0, 2ir] define a point on the 
unit three-dimensional Bloch sphere. The expression of 
Qo w as obtained in Ref. [l6j], fixing the bias of the coin 
toss 9 — 7r/4, following the method developed by Nayak 
and Vishwanath [29j 

Qo = — (1 p) cos 7 + sin 7 (cos p + 1V2 sin 9?) 

2 v2 L 

(32) 

Using this result in Eq. (fT5)l the dependence of x with the 
initial conditions is given by 
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FIG. 2: Isothermal curves as functions of the dimension- 
less angles 7 and ip. Due to the rotation symmetry in 
the angle ip, only four zones are distinguished: two "cold" 
and two "hot". The "hot" zones (orange on line) have six 
isotherms; from inside to outside their temperatures are: 
T/Tq = 6.5, 3.2, 2.2, 1.6, 1.3 and 1.1. The "cold" zones 
(blue on line) have five isotherms; from outside to inside their 
temperatures are: T/T = 0.9, 0.8, 0.7, 0.68 and 0.66. The 
straight (green) lines corresponded to T/Tq — 1, see Eq. H35p . 



where xo = 3/4- 1/V2. 

It is useful to define a characteristic temperature (in 
units of e) 



To = 2/ 



In 



1 



(34) 



X = Xo (1 + cos ip sin 27) , 



(33) 



in order to express any other temperature as a propor- 
tionality with To = 1/Po- Then from Eq.((27|) we obtain 
an expression for j3 as a function of the angles 7 and ip 

( tanh B \ 2 
coswsiii27 = — — - — 1. (35) 

^ V taim /W 

Figures [5] and [3] show the level curves (isotherms) for the 
entanglement temperature as a function of the QW initial 
position. In Fig. [5] the initial position is defined through 
the angles 7 and ip and in Fig. [3] it is defined through 
the position on the Bloch sphere (see Eq.([3"Tj)). Both 
figures show four regions, two of them corresponding to 
temperatures T > T (orange color on line) and the other 
two to temperatures T < T (blue color on line). The 
longest isotherms (green color on line) correspond to the 
temperature T = T n and their initial conditions are 7 = 
0, tt/2, 7T and (p — tt/2, 3tt/2. 
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FIG. 3: The isothermal curves of Fig. [2] shown on the Bloch 
sphere. The QW initial chirality determines the entanglement 
temperature, see Eqs. (|27l31l33fl . 



V. DISTRIBUTED INITIAL CONDITIONS 

In previous works @, H(| we have studied the QW with 
extended initial conditions. Now the entanglement tem- 
perature is studied in such a case. The following extended 
Gaussian distributions is proposed: 
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1 



do V 2vr 



exp 



1 



<7ov 2ir 



exp 



2al 



2a\ 



cos( 7 /2), (36) 



sin( 7 /2), (37) 



where oq is the initial standard deviation, 7 € [0,7r] de- 
termines the initial proportion of the left and right chiral- 
ity and <p £ [0, 2ir) is a global phase. Using these initial 
conditions, Eqs. (j3"6l I3"?)) . the asymptotic value of Q(t), 
see Eas. (|6l7p . was obtained Q as 



1 



with the restrictions 



• cos 7 tan 9 , 



(Jo > 1, 



and 



tan 9 



COS if 



tan 7 

Replacing Eq.® in Eq.([l8]) we obtain 

2 



X 



/ cos 7 \ z 
V 2 cos e) ' 



(38) 



(39) 



(40) 



(41) 




FIG. 4: Isothermal curves as functions of the dimension- 
less angles 7 (initial conditions) and 8 (bias of the coin). 
Four curves (with different colors on line) are presented; each 
curve has two branches placed symmetrically on both sides 
of 7 = 7r/2, where T = 00. The values of T are given by 
Eqs. (|42I) . From left to right the values of T are, in units of 
e: 0.5 (purple), 1. (blue), 2. (green) and 5. (orange). The 
diagram has two discontinuities in 9 — tt/2 and in 7 = n/2 
(see Eqs. (|42l43l) b 



and then using Eq. ([27|) we have 



1 j /|cos#| + |cos7| 

2 n V |cos 6>| - |cos7| 



(42) 



where taking into account Eqs.(|20 ll4T|) the initial condi- 
tion satisfies the constraint 



|cos7| < |cos( 



(43) 



The functions Qq, \ and j3 vanish for 7 = tt/2 ( see 
Eqs. ([3^1l4ni4^|) ) and simultaneously the entanglement 
entropy Eq. (|14[) has its maximum value Sq = 1. This 
maximum value is achieved when the entanglement tem- 
perature is T — 00. Under these conditions the system 
behaves as a classical Markov process [7]. On the other 
hand, the initial conditions 7 and ip are not independent 
(see Eq. (|40"|) ) and for each value of 7 there is only one 
value of T, then for fixed 9, it is not possible to have 
isotherms as functions of 7 and ip. Instead the entangle- 
ment temperature depends on 9 and 7 from Eq. (|42[) , Le- 
the choice of the bias of the coin toss 9 or of the initial 
proportion of the chirality 7 could lead to the same en- 
tanglement temperature. Fig. U shows the isotherms as 
functions of 7 and 9. 
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VI. TRANSIENT BEHAVIOR 
0.08 -i 



0.04- 




t 



FIG. 5: Envelope of the probability A+ — A+ as a function of 
the dimensionless time t for two different initial conditions. 
Each initial condition is established by the couple (ip, 7) in 
Eq. (|3"T)l . Their values are (71-/8, 7r/4) for the full black line 
and (7r/4, 7r/3) for the dashed red line. In both cases, the 
temperature is T = 0.79 lo- 
in the QW a stationary entanglement is established 
between the chirality and position degrees of freedom af- 
ter a transient time. This fact allowed us to introduce 
the concept of entanglement temperature. The transient 
behavior of the system is studied using the original map 
Eq.([3} in a numerical code with initial conditions given 
by Eg. pip . These numerical calculations are summa- 
rized in Figs. [5] and [6] Fig. [5] presents the difference be- 
tween the transient (A+) and the stationary (A+) eigen- 
values of the density matrix as a function of time (see 
Eas.(|13 | ll5 ]) ). The figure only presents the envelope of 
the curves because the real eigenvalues dynamics is very 
intricate; it presents quick oscillations with high density 
of paths where it is only possible distinguish its global 
contour. However, the average evolution of the system 
is determined by the envelope dynamics. Each envelope 
has two branches placed symmetrically on both sides of 
A + — A + = 0. Two pair of curves are presented in full 
and dashed lines with color black and red respectively. In 
both cases, the envelopes decay for t — > 00 as a power law 
l/t c with c = 0.490 for the dashed red line and c = 0.486 
for the full black line. The envelopes of A± will be called 
A± respectively. 

It was numerically verified for several initial conditions 
given by Eq. ([51]) that the transient behavior of A± — A± 
can be adjusted by a power law function of time. Fig. [5] 
shows the power law exponent c as a function of the initial 
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0.46- 



0.44- 



0.2 0.4 0.6 0.8 1 

y/n 

FIG. 6: The power law exponent as a function of the dimen- 
sionless angle 7. The initial conditions (<p, 7) correspond to 
the isotherms T = 1.1 To. 

condition 7 for the same temperature. Remember that, 
for T and 7 given, tp is determined by Eq. (|35[) . There- 
fore, the exponent c has a dependence with the initial 
conditions, however this dependence is not determined 
by the asymptotic temperature value. 

Additionally, the transient behavior of A± — A± was 
numerically studied, using initial conditions given by 
Eqs. (|36H37l) with a >■ 1. In these cases, the system 
showed a negligible transient dynamics, in concordance 
with the calculation developed in Ref. 0- Therefore, 
for these initial conditions, the reduced density matrix is 
essentially always in thermodynamic equilibrium. 

With the aim to understand the transient behavior of 
the system we develop an analytic theory implementing a 
parallelism between the reduced density operator Eq. (|25[) 
and the density operator of an electron in a external mag- 
netic field. With this picture on mind we propose the 
following master equation for the probabilities A+ and 
A_ 

dX - i ^ 
— — = A+ - A_ 

= X-W-+ -A+w + _, (44) 

where W-\ and w |_ are transition probabilities per unit 

of time, that can be understood as population rates, 

corresponds to the transition A+ >— » A_ and w |_ cor- 
responds to the transition A_ >— > A+. These rates are 
time-dependent functions, and their behaviors are known 
in the limit t — > 00 when d\±/dt — > 0. In this limit, 
the stationary solution of Eq. (|4"4"]) must be the couple A_ 
and A_|_, given by Eq. (|15l) . Then the asymptotic values 
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of the population rates satisfy 



VII. CONCLUSION 



w b _ A_ 

W a A + ' 

where w a and Wb are defined by 



Wb 



lim w \- 

t — > oo 

lim W-\ 

t — > oo 



(45) 

(46) 
(47) 



Eq.(|45j) expresses a condition of detailed balance which 
says that the rate of occurrence for any transition equals 
the rate for the inverse transition. Using our knowledge 
about the transient and asymptotic behaviors, the fol- 
lowing population rates are proposed 



= Wb + £(*), 
w-+ = w a - £(t), 



(48) 
(49) 



where 



K c 
= — [w sin(wt + 5) + ( — w a — w b ) cos(ujt + 5)], (50) 
t c t 

with c > 0, K, oj and S constants. The general solution 
of Eq. (|4"4")) with these population rates is 

A+ = A+ + — cos(ujt + 5) + de-^ Wa+w ^\ (51) 
t c 

\_=A_-— C os(ut + 5)- de-( w « +w ^\ (52) 
t c 



where d is 

-{w a +w b )t 



an additional constant. Note that 
e -^ a -r^ o; . _^ o faster than l/t € for t oo, then 
Eqs. (|51l52[) verify the asymptotic behavior obtained nu- 
merically. All the constants K, c, cj, 8 and d depend 
on the initial conditions, however their values should be 
compatible with the positive character of the functions 
A+, A_, W-i and w |-. In particular, to describe cor- 
rectly the numerical results, K takes a finite value for lo- 
calized initial conditions and K takes negligible value for 
distributed initial conditions. In summary, the Brownian 
motion equation for our reduced density matrix Eq. (f2"2"j) , 
takes the form of a master equation Eq. (|4"4")) . 

Finally it is important to point out that the asymp- 
totic behavior found here is similar to the behavior of 
the simple cellular automaton known as sandpile [32| . 
Such behaviors are characteristic of extended dynamical 
systems with spatial degrees of freedom. They naturally 
evolve to self-organized states with correlations that de- 
cay with a power law. 



The unitary evolution of the QW in a composite 
Hilbert space is studied. In particular the entanglement 
between chirality and position degrees of freedom is in- 
vestigated. After a transient time the system establishes 
a stationary entanglement between the coin and the posi- 
tion that allows to develop a thermodynamic theory. The 
asymptotic reduced density operator is used to introduce 
the entanglement thermodynamic functions in the canon- 
ical equilibrium. These thermodynamic functions charac- 
terize the asymptotic entanglement. It is shown that the 
QW initial condition determines the system's tempera- 
ture, as well as other thermodynamic function. A map for 
the isotherms is analytically built for arbitrary localized 
initial conditions. Additionally, it is shown that choosing 
appropriately the bias of the coin-toss, it is possible to 
obtain a predetermined entanglement temperature. 



The transient dynamics of the reduced density opera- 
tor outside the thermodynamic equilibrium is also stud- 
ied. We show numerically that this transient behavior 
can be adjusted with a power law, whose exponent de- 
pends on the initial conditions. We built a master equa- 
tion to describe this behavior where the population rates 
have a time dependence. The accuracy of the master 
equation solution is numerically verified and it is shown 
that the reduced density has a cellular automaton behav- 



The behavior of the reduced density operator looks dif- 
fusive but it has a dependence with the initial conditions, 
the global evolution of the system is unitary. Then, if an 
observer only had information related with the chirality 
degrees of freedom, it would be very difficult for him to 
recognize the unitary character of the quantum evolu- 
tion. In general, from this simple model we can conclude 
that if the quantum system dynamics is developed in a 
composite Hilbert space, then the behavior of operators 
that only belong to one sub-space could camouflage the 
unitary character of the global evolution. 
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